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Abstract 

We investigate the overlap of the ground state meson potential with sets of mesonic-trial wave functions correspond- 
ing to different gluonic distributions. We probe the transverse structure of the flux tube through the creation of 
non-uniform smearing profiles for the string of glue connecting two color sources in Wilson loop operator. The 
non-uniformly UV-regulated flux-tube operators are found to optimize the overlap with the ground state and display 
interesting features in the ground state overlap. 
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1. Introduction 

A fundamental property of the non-perturbative 
regime of confining pure-gauge theories is the linear in- 
crease in the ground state potential between a pair of 
static color sources. In addition to that, lattice gauge 
simulations have recently confirmed the existence of a 
sub-leading non-Coulombic long range correction to the 
mesonic ground-state potential in Yang-Mills theory 
0]. These features are connected with the underlying 
gluonic picture and the subsequent energy distribution 
profile. Although the properties of the ground state po- 
tential have been unambiguously measured to a sub- 
leading order in the infrared region of the non-abelian 
gauge, the geometrical aspects of the associated energy 
density profile at low temperatures remains to be com- 
pletely resolved. 

Lattice calculations of the gluon-field distribution in 
static mesons using Wilson loop operators reveal uni- 
form energy and action density profiles along the line 
joining the static qq pair at large distances [5]. These 
measurements, however, may be vulnerable to system- 
atic errors associated with excited-state contamination 
0]. The non-ground state components manifest them- 
selves in the revealed gluonic profiles as a bias reflect- 
ing the form of spatial links of the Wilson loop operator. 
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The bias by the geometry of spatial links in the L shape 
baryon operator provides a clear example where the flux 
distribution mimics the source latZIl- The height and 
the width of the distribution also depend on the ultra- 
violet properties of the gauge links in the source |@|. 
Apart from the arbitrariness in adopting the source that 
best approximates the ground state, the statistical fluc- 
tuations impose a practical constraint on the Euclidean 
time evolution in the loop operator to isolate the phys- 
ically interesting energy-density profile of the ground 
state. The excited-state contamination is more challeng- 
ing in the case of field distribution calculations which 
involve three-point correlations rather than the potential 
which is extracted in the large time limit of a two point 
correlation 06J] . 

In the finite-temperature regime, the static meson can 
be constructed using a pair of Polyakov lines. These 
hadronic operators provide a systematically unbiased 
stringless gauge-invariant objects in the calculations of 
field-distribution correlations. This means that one need 
not adopt specific geometric or UV properties for the 
gluonic string between the color sources. 

Current investigations of the flux-tube profile in the 
finite-temperature regime of QCD have revealed action- 
densities of non-uniform distribution along the flux-tube 
IH]. The action density displays a two dimensional 
Gaussian-like profile and isosurfaces of a curved pro- 
late spheroid-like shape [8] in the intermediate source 
separation distance region 0.6 < R < 1 fm. This has 
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been observed near the deconfinement point T ~ 0.9T C 
and remain manifest at the temperature T ~ 0.8 T c close 
to the end of the plateau region of the QCD phase dia- 
gram [9]. The measurements of the tube's mean square 
width profile indicate, however, almost constant width 
topology. Variation in the amplitude give rise to curved 
isosurfaces. At larger distances, the tube changes width 
along the qq plane and this width profile is predictable 
based on a free bosonic string picture H. The glu- 
onic distributions obtained at finite temperature by cor- 
relating two Polyakov lines constitute an interesting 
source of knowledge for investigating the possibility 
that non-uniform densities provide the true geometry of 
the ground state in the static meson. The viability of 
considering finite temperature results as an indication 
for the field distribution of the system's ground state can 
be justified by a rguing a lso that the change in the string 
tension is small IMUM at T * 0.87V 

In addition to this observation, it has been found re- 
cently that a model of Coulombic trial states provides 
a good overlap with the ground state in the continuum 
limit 1 12]. Moreover, the free bosonic string model pre- 
dicts observable edge-effects at zero temperature for the 
width profile of the tube given by 113 ] 



— log|cos(^)|,£e [-R/2,R/2] 
Tier R 



(1) 



The above term describes the geometrical shape of the 
flux tube and it indicates subtle changes in the tube's 
mean square radius in the middle of the tube and more 
pronounced changes near the quark positions. The suc- 
cess of the string picture in accounting for the flux-tube 
curvature at high temperature at large distances is re- 
markable, and one may investigate such effects at zero 
temperature. Apart from the string's width effects, a 
non-uniform action density amplitude pattern along the 
tube has been observed at finite temperature whether the 
tube exhibits a non constant width profile or not fl. 
The Bag model is also anther scheme that predicts an 
ellipsoidal-like [ 14] shape for the tube in the infrared 
region. 

At zero temperature, correlating a pair of Polyakov 
lines with an action density operator is very noisy and 
requires substantial numerical simulations using special 
techniques such as the Multi-level algorithm [1]. Never- 
theless, standard Wilson loop operators do not exhaust 
the possibility of investigating the transverse structure 
of the field distribution. We can introduce the non- 
uniformity by employing the idea of constructing the 
flux-tube operator as a product of locally smeared links 
with varying smearing extents. This corresponds to im- 
posing a local transverse cutoff on the lattice parallel 



transporters between the fermionic fields. By extending 
the space of mesonic states constructed in Wilson loop 
this way, we investigate the existence of states such that 
the overlap with the ground state is maximized. This is 
the objective of the present report. 

2. Wilson loop operator 

In the mesonic Wilson loop, a mesonic state is de- 
scribed by fermionic fields connected by the parallel 
transporters Q that is an element of the corresponding 
local gauge group G 



(2) 



The spectral expansion of the Wilson loop operator 
reads 

{W(R,t)) = Y J C n {R)e- v " {R) '. (3) 



The overlaps, C„(R), obey the normalization condition 

(4) 



£c„(fl) = l. 



For large f, the so-called overlap with the system ground 
state can be measured as 



C (R) = 



(W(R,t)Y 



(5) 



(W(R,t + 1))' 

The mesonic state with an infinitesimally thin flux tube 
operator between quarks 



=9 exp 



[ J ' dz. A(z)] 



(6) 



is poorly overlapping with the ground state in the con- 
tinuum 1I12J]. Analogy with an Abelian analytically solv- 
able case B15I1 shows that the infinitesimal thickness of 
the gauge links corresponds to the removal of UV cutoff 
on the transverse direction of the tube causing a vanish- 
ing overlap with the ground state. 
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Figure 1 : Schematic representation of link-blocking 



Applying a smearing operator on a gauge link in a cer- 
tain local gauge group would alter its UV properties. 
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Smearing introduces a transverse UV regulator into the 
flux-tube operator and this results in the enhancement of 
the overlap with the ground state. Define a smearing op- 
erator >S : (B X G — > ® X G. A standard APE II 1611 smearing 
sweep (SU)^(x) consists of a replacement of the spatial 
link-variable U M {x) (yu = 1 , 2, 3) with the angular part of 

U^(x) =(1 - a) U^x) + j J^vWI/pOt + v) 

U'l(x +fi) + U'l(x - P)U M (x - P)U„(x + v - //)}, 

(7) 

where a is the smearing parameter as illustrated Fig. [1] 
In QCD, this corresponds to the projected link U^x) e 
SU(3) C that maximizes 



RTt{&„(jc)^(jc)}. 



(8) 



The geometrical characteristics of smearing can be de- 
scribed by analogy to the Brownian motion associated 
with diffusing a scalar field 01711 . Given a scalar field 
<p(r;n + 1) similar to the (n + l)-th smeared gauge link 
in the fi direction, and a smearing time r = na T with 
a spacing a T . The smearing operation will then corre- 
spond to the diffusion initial value problem, 

5 T 0(r;r)= J D(^+^) < *(r;T), 
0(r; n = 0) =S(r) 

with the diffuseness, 



along the spatial links in Wilson loop. This sequence of 
numbers maps into the geometrical space of the corre- 
sponding radii r(x,) and the amplitudes A(x,) given by 
Eqs. ( [Tol l and (fTTT i. For a mesonic state «,■ or r(x,), the 
projection on the system's ground state is measured as, 



(13) 



The mesonic state constructed by operators correspond- 
ing to a rectangular shape given by a constant sequence 
{r,} has been considered to provide a good approxi- 
mation for the potential ground state ||5|]. The under- 
standing of the geometry of the flux-tube, nevertheless, 
would be increased by constructing trial states with- 
out this constraint. We expect based on the results in 
Refs. S 1211 that the best possible approximation of the 
ground state may be approached this way. 
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Figure 2: Schematic diagram of the smearing profile, 
h is the minimum number of smearing sweeps applied 
at the last link giving rise to smearing radius L\, n is 
the maximum number of smearing sweeps in the middle 
resulting in the radius L2. 
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The Green kernel of the above heat equation partial dif- 
ferential equation, Eq. ©, gives the evolution of the 
scalar field in the smearing time, 



G(r;r) = 



1 



(4ttDt) 



exp 



r.r 



(10) 



The diffuse field is Gaussian distributed with a charac- 
teristic radius, 



'/ d 3 rG(r;i-)r 
/ dhG(r;T) t 



2a^an. (11) 



Applying the local smearing operator S at each spatial 
link in Wilson loop operator, 



= (.S ni U)(x 1 )(S n2 U)(x 1 +d)- 



(12) 



The sequence of numbers of smearing sweeps applied at 
each link {«i, «2, 13, «;v} fixes the gluonic distribution 



For a Wilson loop with spatial links, we consider 
a class of distributions characterized by the radii r = 
f(x) + L\, with an ellipsoidal constraint 



a 2 



b 2 



(14) 



The abscissa x are lattice coordinates and are mea- 
sured from the middle plane between the two quarks 
x = 0. The shape is fixed by the minimum number 
of smearing sweeps at the last spatial link = h and 
the maximum number of sweeps at the middle links 
«A72 = n. The shape consists of a base defined by the 
family of rectangles of height L2 °c yh and ellipsoidal 
caps with a = R/2 and b 2 oc (« - h), thus, the radius 
at middle link is L2 = b + L\. This parametrizes the 
geometrical shape schematically represented in Fig. [2] 
Among a variety of heuristic shapes, this particular pre- 
scription is found to be especially useful for maximiz- 
ing the overlap with the ground state by variations of the 
tuning parameters h and n. 
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heat bath algorithm 11811 using Fabricius-Haan and 
Kennedy-Pendelton (FHKP) 03, updating. Each 



3. Numerical results and discussions 

We take our measurements on 200 SU(3) pure-gauge 
configurations. The configurations are generated using 
the standard Wilson gauge action S w on a lattices of a 
spatial volume of 36 3 x 32 for the considered coupling 
value of/3 = 6.00. 

The lattice spacing at this value is a = 0.1 fm. The 
Monte Carlo updates are implemented with a pseudo- 
sing Fab 

EE EH 

update step consists of one heatbath and 5 micro- 
canonical steps. The measurements are taken on con- 
figurations separated with 1000 updating sweeps. 

The APE smearing operation Eq. (0 and Eq. ((8} is 
locally applied on spatial links of the Wilson loop with 
smearing parameter a = 0.7. That is, the number of 
smearing sweeps at each link as one moves from the 
quark to the antiquark is not necessarily equal. The 
smeared links are drawn from sets of smeared configu- 
rations corresponding to 1 to 40 sweeps of APE smear- 
ing. The spatial links in the Wilson loop are drawn from 
theses sets. 

For noise reduction, the Wilson loop is calculated at 
each node of the lattice and then averaged over the 4- 
volume of the hypertoroid. The overlap with the ground 
state Co of Eq. (f5]lis measured using Wilson loops of 
temporal extent of 2 and 3 slices for source separations 
R = 10a and/? = 12a. 

The sequence of the numbers of smearing sweeps ap- 
plied at each link labels a trial state. Here, we consider 
measurements of Co for states in the parameter space 
{5 < n < 40, 1 < h < 30). The state is uniquely deter- 
mined by n and h. The number of sweeps at each link 
in between is obtained from Eqs. (TT4l and (fTTT i. The 
smearing profile is symmetric with respect to the mid- 
dle point between the quarks. 

Figure[3]indicates the measurements of the overlap of 
the ground state for three selected lines in the parametric 
space correspond to \n,h = 1), \n,h = 5), {n,h = 13). 
The shape is elliptic with the quark source at the end 
of the ellipse for h — 1 (one smearing sweep at the last 
link). The values of Co are small for ellipses with num- 
ber of sweeps in the middle 5 < n < 17. With further 
increase in the height «, the value of Co increases and 
lies approximately in the range [0.80,0.85] for n > 17. 
However, if the height of the rectangular shape in the 
base increases to h — 5 in terms of sweeps, the overlap 
with the ground state is in the range 0.90 < Co < 0.93 
for n — 19 to n = 33. Further increase of h causes subtle 
increases in the value of Co until an optimum value for 
h =i 13 is reached. 



The flux tube modeled as an ellipsoid with the color 
source lies at the far end of the prolate shaped glue does 
not provide the optimal overlap with the ground state 
in static mesons. The subsequent increase in the ground 
state overlap value Co observed by introducing this rect- 
angular base of height h can be understood as consider- 
ing a model of an elliptic-like shaped flux tube in which 
the quark positions are shifted from the edges to the in- 
side. The bag model, for example, leads to an ellipsoidal 
approximation for the gluonic field distribution around 
the color source J3- The above result would indicate 
that the ground state gluonic bag could have the position 
of quarks not exactly at the edges. Indeed the elliptical 
shapes revealed at small quark separations in Ref. 101 
contained the quarks. 

The retrieved values of the overlap for parameter val- 
ues corresponding to the two lines \n,h = 15, 18} are 
illustrated in Fig. |4] The overlap with the ground state 
exhibits a pronounced oscillatory behavior versus n for 
h > 13. Nevertheless, the measured data are not ran- 
domly scattered in the graph. 

The data are seen to arrange themselves to lie ulti- 
mately within what resembles a band structure. This 
branching is more evident when plotting a denser region 
of the parametric space as in Fig. [5] for sweeps n > 25. 
The data appear to line up into five bands with the con- 
tinuous variation of n indicating that this observed oscil- 
latory behavior by changing n for a given h as in Fig. [3] 
may not be arbitrary. This likely to arise from the dis- 
crete nature of «, and h in constructing the source and 
the sink and the inclusion of spatial link configuration 
which systematically probe excited states of the glue. 
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Figure 3: The overlap with the ground state Co, the dis- 
tance between the quark anti-quark source is R — 10 a, 
(3 = 6. 
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There is, however, a variety of states of interest that 
maximize the overlap value at Co ~ 0.94. These states 
line up in the first band from above as in Fig. The 
form of the corresponding operators for four of these 
parametric states are shown in Fig. |6] 




Figure 4: The overlap with the ground state Co, for 
R - 1 fm. The line connects the states corresponding to 
variation of the ellipse semi-major axis for each rectan- 
gular base corresponding to sweeps h — 15 and h = 18. 




Figure 5: Comparison between non-uniformly smeared 
profiles n + h and uniform states n — h represented 
by the smooth line. The smooth line connects central 
values of Co for n = h. The quark source separation 
distance R — 1 fm. 



It is evident that the states with very large values 
of h, for example, (h=24, n=32) tend to assume a more 



flat shape rather than the clear difference in amplitude 
along the tube as in the state (h=13, n=34). Neverthe- 
less, the four operators overlap with the ground state 
equivalently. 

For comparison, the values of Co corresponding to the 
uniformly smeared (flat) states, n = h, is also depicted 
in Fig. [5] In this case, Co is a smoothly varying function 
of n. In addition, the curve interestingly crosses through 
the states of the second band from above. Inspection of 
Fig. [5] shows that states constructed by non-uniformly 
smeared links can maximize the overlap with the ground 
state in a comparative way to the mesonic states with 
uniform flux tubes {« = h). 
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Figure 6: The flux tube operator Eq. (fTOt . Each oper- 
ator consists of family of five Gaussian. The operators 
correspond to the states (h=13, n=34), (h=15, n=32), 
(h=20, n=34) and (h=24, n=33). Theses states maxi- 
mizes the value of the with the ground state. The source 
separation distance R = 10. 



In general, we do observe that neither smearing ap- 
proach can overlap higher with the ground state for all 
the variables of the considered parametric spaces. In the 
far region (large values of n or h), however, there exists 
many states belonging to the highest band in Fig.|5]for 
which the overlap with the ground state, Co, is higher 
than the corresponding flat smearing. This observation 
indicates that smearing near the quark anti-quark pair, 
for sweeps greater than 30, increases the excited-state 
contamination. The links near the quark positions ex- 
hibit different UV behavior from the links at the mid- 
dle with respect to the ground state overlap. The mea- 
sured data provide an explanation why the usual flat 
APE smearing decreases the overlap with the ground 
state for large number of smearing sweeps. 



Recalling that the local smearing operator S"U not 
only alters the width of the fat link at each locus which 
is proportional to yfn, but also decreases the amplitude 
as indicated in Eq. ( TTOb . We see that some of these 
states that maximize the ground state overlap show in 
addition to the variation in the amplitude a sensible 
variation in the amplitude of the flux tube operator near 
the quark source. 
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Figure 7: Same as Fig|5] for quark-antiquark separation 
distance R = 1 .2 fm. 




Figure 8: The overlap with the ground state Co versus 
the physical width of the flux-tube operator in the mid- 
dle Lq of the tube. Each line corresponds to a fixed ratio 
between the width of the tube in the middle and at the 
last link L\, The quark-antiquark separation distance of 
R — 1 .2 fm is considered here. 



At larger source separation R=1.2 fm, the collec- 
tion of Co in branching bands is less obvious as can be 
seen in Fig. [7] The observable difference between the 
overlaps value at R — 1.2 fm in Fig. [3] and for R = 1.0 
fm in Fig. |4] is that the states \h = 24, n > 32), having 
a flatter and more uniform profile, do not optimize the 
overlap with the ground state. 

We replot the points in the parametric space which 
correspond to a fixed ratio between the radii at the end 
and middle points of the flux tube Z4/L2 in Fig. [8] The 
lines pass through trial states that approximately have 
the same non-trivial topology for the case L\ + Li- The 
flat states, L\ — L2, are smoothly varying in comparison 
to the non-uniformly smeared states which in general 
assume the same behavior except for an obvious exis- 
tence of fluctuation along each curve. These arise due 
to the fact that non-uniform states, because of the lat- 
tice structure, do not have exactly the same non-trivial 
topology in comparison to the flat states. 

4. Conclusion 

The overlap with ground state static mesons has been 
measured for a variety of trial mesonic states corre- 
sponding to non-uniform gluonic distributions. An opti- 
mal ground state overlap for non-uniform flux-tube op- 
erators as well as flat smeared operators has been found. 
This supports the possibility that the true ground state 
flux tube is not uniform but rather has a curved flux 
strength profile larger in the middle with higher action- 
density suppression. Such a result resembles the pro- 
file revealed at finite temperature. The findings of this 
work motivate the use of a Multi-level approach to ex- 
plore the action-density profile of a static meson at zero 
temperature. This is the focus of the forthcoming inves- 
tigation. It is remarkable that highly non-uniform trial 
states produce values of Co equally as good as the tradi- 
tional uniform smearing approach. This result indicates 
that it is the smearing extent is the most critical in ob- 
taining optimal overlap with the ground state. However, 
it is important to note that the shape of the non-uniform 
sources is critical to obtaining large values for Co. Our 
investigation of other source shapes complementary to 
Fig.[2]did not produce values for Co as large as favorable 
case studied here in detail. 
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